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To calculate the intensity of x-ray emission 
in electron beam microanalysis requires a 
knowledge of the energy distribution of the 
electrons in the solid, the energy varia- 
tion of the ionization cross section of the 
relevant subshell, the fraction of ioniza- 
tions events producing x rays of interest and 
the absorption coefficient of the x rays 
on the path to the detector. The theoretical 
predictions and experimental data avail- 
able for ionization cross sections are limited 
mainly to K shells of a few elements. 
Results of systematic plane wave Born ap- 
proximation calculations with exchange 
for K, L, and M shell ionization cross sec- 



tions over the range of electron energies 
used in microanalysis are presented. Com- 
parisons are made with experimental 
measurement for selected K shells and it is 
shown that the plane wave theory is not 
appropriate for overvoltages less than 2.5 V. 
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1. Introduction 



The intensity of x rays emitted when an electron 
beam strikes a sample depends on the energy distribu- 
tion of the electrons in the solid, the energy variation of 
the ionization cross section of the relevant subshell, the 
fraction of ionization events that give x rays in the line of 
interest and the absorption coefficient of the x rays on 
the path to the detector. This can be summarized as [1] 



/(£, r)(7x(£)/xexp(-/Axlro - rl)drd£: (1) 



where /(£", r) is the distribution of electrons in the spec- 
imen as a function of energy E and position r, axiE) is 
the ionization cross section for the relevant subshell, fx 
is the fraction of ionization events producing x rays in 
the line of interest and Ox is the absorption coefficient 
for the X rays on their path to the detector at position Tq. 
The integration is over all electrons that can ionize the 
subshell of interest, energy Ei, and over the volume of 



the specimen. For TEM this expression can be consider- 
ably simplified as specimens are so thin that energy loss 
is negligible. (The mean energy loss can actually be 
measured from the energy loss spectrum which could be 
recorded at the same time as the x-ray spectrum in 
suitably equipped microscopes. A typical range of val- 
ues for the mean energy loss is about 20 eV to 50 eV 
which is very small compared to the microscope accel- 
erating voltage.) The cross section is then just the ioniza- 
tion cross sections for electrons at the beam energy. 

In microanalysis the ionization cross section for pro- 
duction of X rays and the electron energy distribution are 
often multiplied together to form a new function <P(pz), 
that is a function of depth in spatially homogeneous 
specimens [1]. 



Ix 



(r)exp(- fjixlro - r\)dr 



(2) 
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(f>{r)= I{E,r)ax{E)f^dE 



(3) 



There has been much debate in the Hterature on the 
shape of the <P(pz) function [2], A reUable standardless 
scheme for determining elemental composition from x- 
ray intensities requires that all these processes be mod- 
eled accurately from the relevant physical theory or val- 
ues tabulated from experimental measurements. 

The energy distribution of electrons in the specimen 
can be modeled either from numerical solutions of the 
Boltzmann transport equations [3] or Monte Carlo cal- 
culations [4]. Each method has its strengths and weak- 
nesses. The methods that Fathers and Rez [3] used for 
numerical solution of the Boltzmann transport equation 
are fast and efficient, but are limited in the number of 
energy levels that can be used. They can be generalized 
to layered specimens but are inapplicable for general 
inhomogeneities. Monte Carlo calculations can be used 
with arbitrary specimen geometries, but care should be 
taken that the sampling is done correctly and that simple 
approximations such as continuous slowing down do not 
lead to significant error. In other fields such as Medical 
Physics a combination of small angle transport theory 
and Monte Carlo calculations is used [5]. 

The fraction of events giving x rays in the line of 
interest is a product of the fluorescence yield and the 
fraction of x rays in a particular line. Reliable calcula- 
tions for fluorescence yield and x-ray emission rates 
based on atomic physics are readily available [6,7,8,9]. 
Absorption coefficients, which could be calculated from 
first principles expressions for the photoelectric effect, 
have also been tabulated. Arguably the least known in- 
gredient in the expression for x ray intensity is the elec- 
tron ionization cross section. Scofield [10] has tabulated 
K and L shell ionization cross sections for a selection of 
elements at relativistic energies that are above the ener- 
gies used in microanalysis. Rez [11] published a number 
of K shell ionisation cross sections based on hydrogenic 
wavefunctions, and L and M cross sections based on 
Hartree-Slater wavefunctions. There have been some 
calculations for K shells for inert gases and selected 
transition elements. 

In practice microanalysts use cross sections based on 
the Powell [12] parameterization of the Bethe formula 



, 2, 6.51 X IQ-'Vfcx , 
crx(cm ) = ^TT. logc 



Cx£( 



(4) 



where n^ is the number of electrons in the subshell and 
bx and c^ are parameters to be determined by fitting to 
experimental data or other calculation. An obvious 
shortcoming with this expression is that it fails for elec- 
tron energies less than CxEi when c^ is greater than 1 



(which is usually the case). The same problem arises in 
the expression for the stopping power and Joy and Luo 
[13] have proposed a modification for these low ener- 
gies. My early work [11] was aimed at using first prin- 
ciples ionization calculations to estimate the parameters 
in the Bethe-Powell expression. 

I have now completed systematic calculations of the 
electron ionization cross sections from all subshells that 
might be relevant for microanalysis. The plane wave 
Born approximation with the Ockhur [14] approxima- 
tion for exchange was used. The lower binding energy 
limit was determined from the lowest energy x-ray line 
that could be detected using an ultra thin window or 
window less detector, that is about 150 eV. In practice 
this means L3 ionisation cross sections were tabulated 
from sulfur, and M5 cross sections from strontium. The 
high binding energy limit was set at about 40 kV which 
meant that K cross sections were calculated up to 
promethium and all L3 and M5 cross sections up to 
uranium were tabulated. The cross sections were calcu- 
lated for a range of energies from the binding energy or 
5 kV, whichever was lower, up to 400 kV. This repre- 
sents the range of electron energies used in both scan- 
ning electron microscopes and transmission electron 
microscopes. The steps were chosen to correspond to 
typical step sizes on electron microscopes, though of 
course intermediate values could be estimated by inter- 
polation. 

A selection of the results is presented and they are 
compared with experimental measurement and other 
previously reported calculations. 



2. Theory 

The complete theory is given elsewhere [15] so in this 
paper only the essential points will be summarized. An 
electron with energy E ionizes an inner shell with bind- 
ing energy Eb, and in the process is scattered by a wave 
vector q emerging with energy E' and ejecting an elec- 
tron of energy s from the atom. Conservation of energy 
requires that 



E' = Eu + e 



(5) 



Since electrons are indistinguishable there is no way to 
make a distinction between the scattered electron and 
the ejected electron so exchange has to be explicitly 
taken into account. The differential ionization cross sec- 
tion is given by [16] 



dV 
d/}dfi 



2y' 
flo 



\f+g\ 



^ + 71/- 



g\' 



(6) 
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where 

/= n^'f(ri)v>'f(r2)7 



1 



Iri - rzl 

is the "direct" term and 

J 

I/-1 - r^l 



<p,(r,)(A,(r,)dV,dV, (7) 



g = nff(r2)<p'f(ri} 



(p,(r,).A.('-.)dVidV, (8) 



is the "exchange" term where the ejected electron ap- 
pears to have changed places with the scattered electron. 
The wavefunction i/'i(ri) represents the incident electron, 
q>i(r2) the inner shell electron, t///(ri) the scattered elec- 
tron and (Pfir2) the ejected electron. The constant ao, the 
Bohr radius, gives the scale of the interaction and -y is the 
relativistic correction factor. Expanding Eq. (6) gives 



dV 
dfids'' 






.f + ; 



1 



k 



(9) 



where the first term in the square brackets is the direct 
scattering which is much larger than the other two 
terms. 

The incident wave i/'i(ri) and the scattered wave ^j{r\) 
can be represented by plane waves exp(/ft:ir) and 
exp(/A:fr) respectively. The momentum transfer, q, is 
just the difference between the initial and scattered 
wavevectors. 



q = k. 



(10) 



with limits defined by the kinematics of the scattering 
[17] 



qmm = Vfei^ + kf - Ikjif 



^n,„ = \/k^ + kf + 2kikf 



where k, the scattered wave vector is 



(11a) 
(lib) 



kf 



nioc 



\E - £b - s) 


2 

+ 2 


'(E - £b - e) 



(lie) 



For simplicity just consider the direct term which can be 
written as 



dV _ Ay'- kf l/y/ (r')exp(/g • r ')yi(r')dVP 

.2 ,, ^4 u^; 



dflAe al h 



It is more convenient to express the differential cross 
section in terms of a quantity known as the Generalized 
Oscillator Strength, F{q, e) defined by 



F{q, ^) = 1^ i;(p;(r')exp(/? • r')(p,(r')dVf (13a) 



such that 



dV 
d/2dfi 



4y' k, h^ 
al k\ 2msq 



2F(q,s). 



(13b) 



To calculate the ionization cross section for the subshell 
Eq. (13b) has to be integrated over all the allowed 
wavevectors and over all possible energies of the ejected 
electron. In our calculations inner shell wave functions 
were taken from the Dirac-Slater program of Liberman 
et al. [18], and continuum wavefunctions for the ejected 
electron were calculated using the self-consistent atomic 
potential. This potentially laborious computation can be 
minimized by examining the behavior of the GOS as a 
function of q and e, known as a Bethe surface [17]. In 
the limit where the ejected electron energy s is about 4 
Eb the ejected electron can be represented by a plane 
wave with wave vector q . The GOS tends to a Gaussian 
shape known as the Bethe Ridge with a half width in 
nm ' approximately given by 35V £3 where Eb the 
binding energy is in Hartrees. The GOS values were 
tabulated to a momentum transfer of 110V£b nm"' 
over an energy range of 4 Eb- Above this energy the 
scaled form of the Bethe Ridge was used. It was found 
that exponential grids in both wave vector and electron 
energy gave the minimum error in the numerical inte- 
gration, and 64 points in q and 64 for s were sufficient 
to assure numerical convergence to better than 2 %. 

An approximate expression for the exchange scatter- 
ing due to Ochkur [14] also assumes that the ejected 
electron is a plane wave. The exchange contributions to 
the scattering can then be formulated in terms of the 
GOS at the Bethe Ridge, ^3(9, e). The GOS in Eq. 
(13b) is replaced by 

q'^ 1 q^ 

F'=(q,s) = F(q, s) + -0FB{q,s) ---^FB(q, e). 

(14) 

This expression only applies when the momentum trans- 
fer is very much less than the scattered electron 
wavevector, which unfortunately does not apply at low 
energies above threshold. Note also that exchange al- 
ways reduces the GOS and hence the cross section. 



3. Results and Discussion 

It is well known that the ionization cross section typ- 
ically rises from threshold to a maximum at about 3 Eb 
and then slowly falls off [12,17]. Since this behavior is 
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universal it is often convenient for comparison purposes 
to plot the cross section against the overvoltage U , the 
ratio of the electron energy to the binding energy. 



U = - 



(15) 



In Fig. 1 the cross sections for potassium K, cadmium L3 
and uranium M5 are plotted as a function of electron 
energy. The dashed curves are the direct contribution 
only, the solid curves incorporate the effects of ex- 
change with the approximate theory. Exchange not sur- 
prisingly has a larger effect near threshold and typically 
can lower the ionization cross section by 20 % in the 
region of the peak at about 3 E-r. At higher energies 
typical of transmission electron microscopes the ex- 
change correction is less important, making a 3 % dif- 
ference for inner shells with binding energy of about 3.5 
kV, and an 8 % difference for inner shells with binding 
energies of about 11.5 kV. The corresponding figures for 
a 400 kV accelerating voltage are 1 % and 3 %, respec- 
tively. 

Since this is the first time that M cross sections have 
been explicitly calculated it is interesting to compare 
them with K and L cross sections for subshells with 
comparable binding energy. The cross sections for 
potassium K cadmium L3 and uranium M5 are all shown 
as Fig. 2, all scaled to the potassium K cross section. 
The shapes are very similar though there are some dif- 
ferences in the peak region. 

There have been very few systematic calculations of 
cross sections reported in the literature, the most com- 
prehensive being those published by Scofield [10]. He 
was mainly interested in extreme relativistic energies so 
there is a very limited overlap with the results presented 
here. Figure 3 shows a comparison for yttrium K . There 
is excellent agreement except at the lowest energies 
where Scofield's cross section is higher. It is hard to see 
whether the discrepancy is due to numerical problems as 
Scofield never shows the interesting region around the 
peak of the cross section. 

While comparison with theory is limited to very high 
energies, the experiments are confined to low energies, 
very often not even covering the peak of the cross sec- 
tion. Nearly all the results are for K shell ionization of 
transition element thin films, no doubt because they are 
easy to prepare. There are additional complications 
from backscattered electrons generating x rays when the 
film is lying on a substrate. Luo et al. [19], have used a 
transport equation treatment to make appropriate cor- 
rections to their measurements. The majority of the 
measurements have been made with dedicated scatter- 
ing chambers, only the recent work of Llovet et al. [20] 
made use of a microprobe. In Fig. 4 the present 
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Fig. 1. lonisation cross sections for (a) potassium K, (b) cadmium Li 
and (c) uranium M5. The solid line with filled circles represents the 
direct contribution, the dashed line with filled diamonds shows the 
result after taking into account exchange. 
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1.510 




Overvoltage 

Fig. 2. Comparison of ionisation CdLs (short dashed line) and UM3 (long dashed 
line) cross sections scaled to the potassium K (solid line) ionisation cross section. 
These subshells all have approximately the same binding energy. 
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Fig. 3. Comparison of the YK calculation with the results of Scofield. The long 
dashed line is without exchange, the solid line is with exchange, the dashed line 
with diamonds are Scofield's results. 



calculations, with and without exchange, are compared 
to various measurements for TiK, CrK, CuK. There is 
very little agreement between the experiments for TiK 
and our calculation lies between the two sets of mea- 



surements [21,22]. For CrK Llovet et al. [20] and Luo 
[19] agree with each other, while He's results [21] are 
much lower. The experimental measurements are in 
broad agreement for copper K, though Llovet et al. [20] 
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Fig. 4. Comparison of experimental results with calculation for Ti, Cr, Cu, and Cu (low energy region). 
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is systematically higher than He [21] and Shima [23]. 
Davis et al. [24] measured cross sections for higher 
energies and these are in good agreement with our cal- 
culations. In the low energy region our calculations for 
both copper K and chromium K underestimate the cross 



section in the low energy region and appear to put the 
peak at too high an energy. 

In his work on the relative cross sections for electron 
and positron emission Hippler [25] argued that the inci- 
dent particle was either accelerated or decelerated by the 
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nuclear electrostatic field. He corrected the incident en- 
ergy by the value of the nuclear electrostatic potential at 
an average radius for the inner shell wave function. Hip- 
pier [25] assumed hydrogenic wavefunctions since he 
was only interested in K shells. His procedure can be 
generalized by calculating the average radial value of 



the wave function explicitly using the wavefunctions cal- 
culated by the Dirac-Slater or Hartree-Slater program 
and then finding the total potential for this point. This 
potential is then subtracted from the incident energy, 
which shifts the calculated curves to the left. The com- 
parison with experiment is plotted again as Fig. 5 for the 
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Fig. 5. Comparison of experimental results with calculation incorporating acceleration of incident electron for 
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low energy region of tinanium K, chromium K and 
copper K and although the peak position is now shifted 
to approximately the correct energy the overall behavior 
in the low energy region is still incorrect. 

The problems arise from an incorrect treatment of the 
wavefunctions for low energies above threshold. If it is 
wrong to treat the ejected electron as a plane wave for 
energies less than 4 ^b, then it is inconsistent to treat the 
scattered wave as a plane wave at these energies since 
electrons are indistinguishable. A correct formulation 
would represent all the electrons as radial solutions of 
the appropriate energy of the Shrodinger equation with 
the atomic self-consistent potential. The theory would 
look very different and the direct them would then be 



X 



/= (4tt)' 



Z3 k I2 
-ni} q OT2 



f) 2 



U k l] \(l4 k l\ 
-OT4 q mj\0 



h k /2\ V(2/4+l)(2/3+l)(2Z2+l)(2/i+l) 
0/ {2k+\) 



x\R,^(s,n)R,^iT- 



, 2r* . 
s, ri) k+iJi. 
' > 



/i k U 



{kr2)Ri,{£, ri)dridr2 
(16) 



where the terms are Wigner 3/ symbols 

\— mi q m2/ 

and r> and r< mean the greater or lesser of ri and r2, 
respectively. The exchange term, g, is identical apart 
from switching ri and r2 in 7?/, (e, rO and R,^ {T — e, r2), 
respectively. This expression is much less amenable to 
fast computation. 



4. Conclusions 

A selection of results from a comprehensive set of 
calculations for K, L, and M shell ionization cross sec- 
tions relevant for microanalysis has been presented. The 
theory is based on the plane wave Born approximation 
with the Ochkur [14] expression for exchange. The 
asymptotic behavior at the Bethe Ridge was used to 
simplify the calculation and to provide values for the 
exchange integral. 

The cross sections have been tabulated over an energy 
range from 5 kV to 400 kV appropriate for both scan- 
ning and transmission electron microscopes. The cross 
sections show similar behavior with energy, showing a 
peak at about 3 E-s- There is good agreement with both 
Scofield's [10] calculations and the small number of 
experimental measurements above this peak. They can 
therefore be used with confidence in transmission elec- 



tron microcopy with beam energies of 100 kV and 
above. The calculations underestimate the value of the 
cross section in the low energy region and show the 
position of the peak at too high an energy. This dis- 
crepancy is due to the use of a plane wave for the 
scattered electron. A product of a radial wavefunction 
and spherical harmonics would be more appropriate for 
the scattered electron wavefunction and a different theo- 
retical formulation should be used in this region. 
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